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Abstract. We evaluate the determinant deti<i^j<„ {^{i^^i_^2j) ~ ix+i+2j))' which gives the 
number of lozenge tilings of a hexagon with cut off corners. A particularly interesting feature 
of this evaluation is that it requires the proof of a certain hypergeometric identity which we 
accomplish by using Gosper's algorithm in a non-automatic fashion. 



The purpose of this paper is to provide a direct evaluation of the determinant 

i<ij<n \ \x-i + 2jj \x + i + 2jJ ) ' 

This determinant arises in our study [4] on the enumeration of lozenge tilings of hexagons 
with cut off corners. For example, consider a hexagon with side lengths x + n, n, y, x + n, 
71, y (in cyclic order) and angles of 120° of which two adjacent corners are cut off as in 
Figure l(a).^ Figure 1(b) shows a lozenge tiling of this region, by which we mean a tiling 
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Typeset by Am^-T^ 



1 



2 



M. CIUCU AND C. KRATTENTHALER 



by unit rhombi with angles of 60° and 120°, referred to as lozenges. The number of these 
lozenge tilings is given by the determinant (1). This is seen by converting the lozenge tilings 
into families (Pi, P2, . . . , Pn) of nonintersecting lattice paths consisting of positive unit steps, 
where the path Pi runs from (z, —i) to (a; + 2z, y — z), z = 1, 2, . . . , n and does not cross the 
diagonal y — x — 1 (see Figure 2), and then applying the main theorem of nonintersecting 
lattice paths [18, Lemma 1], [8], [23, Theorem 1.2] (see [4] for details and background; there 
is also another case in [4] in which the determinant (1) provides the solution). 



X + n 




X + n 

(a) A hexagon with cut off corners. 




(b) A lozenge tiUng of the hexagon with cut off corners. 



Figure 1 



By Theorem 1 below, the number of the lozenge tilings of the preceding paragraph is 
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(a) The path family corresponding to a lozenge tiling. 




(b) The paths made orthogonal. 



Figure 2 

given by a closed form expression. The proof of Theorem 1 that we present in this paper^ is 
primarily based on hypergeometric series identities. A remarkable aspect is that it contains 
an instance of a non- automatic application of Gosper's algorithm [9] (see also [10, §5.7], [20, 
§11.5]), see Step 3 of the proof of Theorem 1. This is noteworthy, because Gosper invented 
his algorithm to automate summation, so that a non-automatic application must be almost 
considered as a misuse. But clearly (and more seriously), the fact that Gosper's algorithm 
is also useful in "computer-free territory" only adds to its value. (The only other instance 



■^An alternative proof is presented in [4], in which a combinatorial argument is used to convert the 
determinant (1) into a different determinant that was already known from [12, Theorem 10]. 
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of a non-automatic application of Gosper's algorithm that we are aware of appears in [19]. 
However, the purpose of use there is different. Roughly speaking, we use it to prove a positive 
result, namely to verify the truth of an identity between certain hypergeometric series, see 
(14). In contrast, Petkovsek and Wilf use it to prove a negative result, namely that a certain 
binomial sum cannot be expressed in terms of closed form expressions.) 

Theorem 1. Let n be a positive integer, and let x and y be nonnegative integers. Then the 
following determinant evaluation holds: 

det ((^+y+j)^(^+y+^ 

^<i,j<n \\x — i + 2j J \x + i + 2j 

^ A (j -l)\{x + y + 2j)\ {x-y + 2j + 1), {x + 2y + 3j + l)n-j , . 
f}^ ix + n + 2j)\{y + n-j)\ ' ^> 

where the shifted factorial {a)k is defined by (a)^ := a{a + 1) ■ ■ ■ {a + k — 1) , k > 1, and 
(a)o := 1. 

Remark. We formulate Theorem 1 only for integral x and y. But in fact, with a generalized 
definition of factorials and binomials (cf. [10, §5.5, (5.96), (5.100)], Theorem 1 would also 
make sense and be true for complex x and y. 

Proof. We prove the determinant evaluation by "identification of factors," a method that is 
also applied successfully in [2], [3], [5], [6], [7], [11], [12], [13], [14], [16], [17] and [21] (see in 
particular the tutorial description in [15, §2.4] or [13, §2]). 

First of all, wc take appropriate factors out of the determinant. To be precise, we take 
{x + y + (^{x + n + 2jy. {y + n — jy.^ out of the j-th column of the determinant in (2), 
J = 1, 2, . . . , n. Thus we obtain 

■Q (x + y + jy. 



fj^ {x + n + 2jy{y + n-jy 
X det {{x + 2j - i + l)n+i{y + i- j + l)n-i - {x + 2j + i + l)n-i{y - i- j + l)n+i) 

l<i,j<n 

for the determinant in (2). Let us denote the determinant in (3) by Dn{x,y). Comparison 
of (2) and (3) yields that (2) will be proved once we are able to establish the determinant 
evaluation 

Dn{x,y) 

= det {{x + 2j -i + l)n+i{y + i-j + l)n-i-{x + 2j + i + l)n-i{y-i-j + l)n+i) 

l<i,j<n 

n 

= ll{j -'^V-ix + y+j + l)i {x-y + 2j + l)j {x + 2y + 3j + l)n-j. (4) 

For the proof of (4) we proceed in several steps. An outline is as follows. In the first step 
we show that nj=i(^ ~ 2/ + 2j + l)j is a factor of Dn{x, y) as a polynomial in x and y. In 
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the second step we show that 11^=1 + y + j + l)j is a factor of Dn{x, y), and in the third 
step we show that 11^=1 + 2y + Sj + l)n-j is a factor of Dn{x,y). Then, in the fourth 
step we determine the maximal degree of Dn{x,y) as a polynomial in x, and the maximal 
degree as a polynomial in y, which turns out to be n(3n + l)/2 in both cases. On the other 
hand, the degree in x, and also in y, of the product on the right hand side of (4), which 
by the first three steps divides Dn{x, y), is exactly n(3n + l)/2. Therefore we are forced to 
conclude that 

n 

Dn{x, y) = C{n) Y[{x-y + 2j + l)j {x + y + j + l)j {x + 2y + 3j + !)„_,■, (5) 

where C{n) is a constant independent of x and y. Finally, in the fifth step, we determine 
the constant C(n), which turns out to equal 11^=10 ~ !)'• Clearly, this would finish the 
proof of (4), and thus of (2), as we already noted. 

Step 1. 11^=1 ~ ?/ + 2j + l)j is a factor of Dn{x,y). Let us concentrate on a typical 
factor (x — y + 2j + I) , 1 < J < n, 1 < / < J. We claim that for each such factor there is 
a linear combination of the columns that vanishes if the factor vanishes. More precisely, we 
claim that for any j, I with 1 < j < n, I < I < j there holds 

I — I 

{j -I) ij + l-2s + {x + 2j + l + n-s + l)s-i {x + n + 2s + l)j+i-2s 

^ (j-s) {s-l)l {2x + 2j + l + s + l)j_s 

■ (column s of Dn{x, x + 2j + /)) 

+ (column j of Dn{x, x + 2j + I)) = 0. (6) 

To avoid confusion, for j — I it is understood by convention that the sum in (6) vanishes. 
In order to verify (6), we have to check 

I — I 

(j -I) U + l- 2s + l)s-i {x + 2j + l + n-s + l)s-i {x + n + 2s + l)j+i-28 
^ {j-s) {s-l)l {2x + 2j + l + s + l)j-s 

■ {{x + i + 2j + 1 - S + l)n-i (x -i + 2s + l)n+i 

-{x-i + 2j + l- s + l)n+i {x + i + 2s+ l)n-i) 

+ {x-i + 2j + {x + i + j + l + l)n-i 

-{x + i + 2j + {x-i+j + l + l)n+i = 0, (7) 

which is (6) restricted to the i-th row. The exceptional case j — I can be treated immediately. 
By assumption, the sum in (7) vanishes for j = I, and, by inspection, also the other two 
expressions in (7) vanish for j = I. So it remains to establish (7) for j > I. In terms of the 
standard hypergeometric notation 
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this means to check 

{x + i + 2j + l)n-i {x -i + 2l + l)n+i+j-l 



{2x + 2j + 21 + 
X4F3 



^ + |, I - i + |, -i - 2i - X, 1 + 2j + 2/ + 2x 



(x -i + 2j + l)n+i {x + i + 2l + l)-i+j-i+n 



(2x + 2j + 21 + 

X 4F3 



2 ^2'2 2^2' 



2j -x,l + 2j + 21 + 2x 



+ I, ^ + ^ + 1 + ^,1+^ + 1 + 



+ (x -i + 2j + l)n+i {x + i + j + 1 + l)n-i 

-(x + i + 2j + ix-i+j + l+ l)n+i = 0. 



(8) 



Both 4F3-series can be summed by means of a 4F3-summation which appears in a paper by 
Andrews and Burge [1, Lemma 1] (see [12, Lemma A3] for a simpler proof), 



4-^3 



N_ 1 
2 ' 2 



N 



„ , a, a + 6 ^ 
^ : 1 



I - ]\[ b lib 
J. lY , 2 ) 2 ^ 2 



(a + 6) TV _^ (-a)jv 



AT 



N 



where A?" is a positive integer. We have to apply the case where N = j — I. This is indeed a 
positive integer because of our assumption j > I. Some simplification then leads to (8). 
This shows that n^=i(^ ~ 2/ + 2j + l)j divides D^ix, y). 

Step 2. 11^=1 + y + J + l)j is a factor of Dn{x, y). Let us concentrate on a typical 
factor [x + y + j + 1), 1 < j < n, 1 < / < j. We claim that for each such factor there is a 
linear combination of the columns that vanishes if the factor vanishes. More precisely, we 
claim that for any j, I with 1 < j < n, 1 < I < j there holds 



E 

s=l+j-l 



J-s 



l-l 

s + l-j-1, 

(x + n + 2s + 1)2,-2. {2x + 3J + 1 + S + 
{x + j + s + |)j-s {x + j + l + s)j-s {x + j + 1 - n + s) 



3-s 



(column s of Dn{x, —x — j — I)) = 0. (9) 



In order to verify (9), we have to check 



E 

s=l+j-l 



4 



l-l 
s+l-j-1 

{x + n + 2s + l)2j-2s i2x + 3j + l + s + 1) 



{x + j + s + I) j-s {x + j + l + s)j-s {x+j + l-n + s)j-s 

■ ((-X +i - j - I - S + l)n-i {X -i + 2S + l)n+i 

- {-X - i- j- l- s + l)n+i {x + i + 2s + l)n-i) = 0, 
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which is (9) restricted to the i-th row. Equivalently, using hypergeometric notation, this 
means to check 



(-1)' 



(- 



2j 



X 



{3 + i + 2j-2l+ x)-2-^+2l+n (2 + 4j + 2x)i-i 



4^-1(1 + 2i + x)i-i (I + 2j - Z + (1 + 2j 



n + xji- 



X 4-F3 
-(-1) 



l,^+2j -l + x,l + 2j + x,l + i + 2j + x 



2 + 4 j + 2a; , 2 + f + j - Z + f , I + I + j - Z + f 
I {i - 2j - x)n-i {S-i + 2j -21 + x)-2+i+2i+n (2 + 4j + 2x)i-i 



4'-i(l + 2j + x)i-i {^+2j-l + x)i-i {l + 2j-n + x)i-i 



X 4-F3 



^ + 2j -l + x,l + 2j + x,l-i + 2j + x,l-l 

3 i \ A 1 I x o i \ A 1 I X 



2 - I + j - Z + f , 2 - f + j - Z + f , 2 + 4 j + 2a; 



= 0. (10) 



In order to estabhsh (10) we apply Bailey's transformation for balanced 4F3-series (see [22, 
(4.3.5.1)]), 



4-^3 



a, 6, c, — 
e,/, l + a + 6 + c — e 

(e - a)N if - o)Ar 



/-iV 



(e)iv(/) 



-4-^3 



AT 



-iV,a, l + a + c- e- /-iV, l + a + 6 



f-N 



l + a + b + c-e-f-N,l + a 



N,l + a- f - N 



where TV is a nonnegative integer, to the second 4F3-series in (10). Thus it is converted into 
the first 4F3-series, and it is routine to check that also the remaining terms that go with 
the 4F3-series agree. So, the two terms on the left hand side of (10) cancel each other, as 
desired. 

This establishes that YTj=i + V + j + divides Dn{x, y). 

Step 3. nr=i('^ + 2y + 3i + l)n-i is a factor of Dn{x,y). This is the most difficult part 
of the proof of (4). Trials of finding linear combinations of columns that vanish resulted in 
extremely messy expressions. So, we decided to work with linear combinations of rows this 
time. Still, the coefficients are not as "nice" as in Steps 1 and 2. 

Let us concentrate on a typical factor {x + 2y + 3i + l), l<i<n, l<l<n — i. We claim 
that for each such factor there is a linear combination of the rows that vanishes if the factor 
vanishes. More precisely, we claim that for any i, I with l<i<n, l<l<n — i there holds 



i+l 



{k + i + l + l)i+i-k 



k=l 



{i + l-k)l 
where Pi (e, /) is the polynomial 



Pi{2i,i + l- k) ■ (row k of Dn{-2y - 3i - l,y)) = 0, 



(11) 



21+1 



^/(e, f)=y^ar ie)r {-f)2i+i- 



r=0 



(12) 



with the expansion coefficients ar given by 



Ur = (x'^) ((a;2 + a; + l)'-^(2x + l)(a; + 2)(a;-l)). 



(13) 
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Here, {x'^)g{x) denotes the coefficient of x"^ in g{x). 

By specializing (11) to the j-th column, splitting the resulting sum into two parts in the 
obvious way, and then moving one sum to the right hand side, we see that in order to verify 
(11), we have to check 

E ^('^ + \ t fc)^"" ^' i + l-k){-2y-3i-l + 2j-k + l)n+k {y + k-j + l)n-k 



k 

i+l 



= E + Z + I + l)i+l-k ^^^^ , + ^ _ fc) (_2^ _ 3, _ / + 2j + + l)„-fe {y-k-3 + l)n+k. 

fc=l 

or, after adding one more term as first summand on both sides, equivalently, 
= E + 2 + / + l),+;-fe p^^^.^ i + i_,,^_2y_si_i + 2j + k+ l)n-k {y-k-j + l)n+k. 



Empirically, we discovered that apparently both sums in (14) are indefinitely summable 
( "Gosper-summable" ; see [10, §5.7], [20, §11.5]). It is exactly this fact which makes (14) 
tractable. 

In the following we will show that the sums in (14) are equal, however, without exhibiting 
an explicit expression for the sums. Instead, what we will do is to read through Gosper's 
algorithm [9] (see also [10, §5.7], [20, §11.5]), which is an algorithm that solves the problem of 
indefinite summation for hypergeometric sums. (For any fixed /, our sums in (14) belong to 
the category of hypergeometric sums.) In the course of reading through Gosper's algorithm 
it will emerge that the sums on both sides of (14) must be equal. 

Let us recall what Gosper's algorithm does and how it works. Let t{k) be a "hyperge- 
ometric term", i.e., be a term such that the ratio t{k + l)/t{k) is a rational function in k. 
Then the Gosper algorithm will find a hypergeometric term T{k) (if it exists) satisfying 



t{k) = T{k + l)-T{k). (15) 

The upshot of this is that then the indefinite summation of the term t{k) can be easily 
carried out, 

B 

Y,m = T{B + l)-T{A). (16) 

k=A 

The term T{k) is found in the following way. First, one finds polynomials p{k), Q{k), and 
r{k) such that 

t{k + l) _ p{k + l) q{k) 
t{k) p{k) r{k + l)' ^ ^ 
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where q{k) and r{k) have the property that whenever (k + a) \ q{k) and {k + (3) \ r{k) 
then the difference a — (3 must not be a positive integer. Next, one finds a polynomial s{k) 
satisfying the recurrence relation 

p{k) ^ q{k)s{k + l) -r{k)s{k) (18) 

for all k. The term T{k) is then given by 

Now let us carry out this program with the summands in (14). First, let t{k) — ti{k), 
where ti{k) is the summand of the sum on the left hand side of (14), 

ti (k) = ^ ^ Pi{2i,i + l-k) (-2y-Si-l + 2j-k+l)^+k {y + k-j + l)n-k. 

Then (17) holds with p(A;) =pi{k),q{k) = qi{k),r{k) = ri(A;), where pi(A;) = Pi{2i,i+l-k), 
qi{k) ^ {i + l- k){-2y - 3i - I + 2j - k), and ri(k) = (i + I + k){y - j + k). So, next we 
have to find a polynomial si{k) satisfying the recurrence 

Pi{2i,i + l-k) = (i + l-k)(-2y-3i-l + 2j -k)si(k+l)-(i + l + k){y-j + k)si(k). (21) 

For each specific instance of i and / this is just routine. However, we were not able to find an 
explicit formula for si{k) in general. Fortunately, we do not need such an explicit expression. 
Assuming that we have found a polynomial si{k) satisfying (21), by (16) and (20) we have 

E i + l-k) i-2y - 3i - I + 2j - k + l)n+k {y + k - j + l)n-k 

fc=o + 

n{i + l + l)si{i + l + l) ri{0)si{0) 

= /■ I , I .S. + 1 + i-) J-T tl(UJ 

pi{i + l + l) pi(0) 

- ^'^^^'+'+'i-2y-3i-l + 2j + l)n{y-j)n+iSi{0), (22) 



the last line being due to the fact that ti{i + I + 1) = 0. 

On the other hand, for t{k) = t2{k), where t2{k) is the summand of the sum on the right 
hand side of (14), 

t2 (k) = ^ \ ^ Pi{2i,i + l-k){-2y-3i-l + 2j + k + 1)^,^ {y-k-j + l)n+k 

we may choose p{k) = P2{k), q{k) — q2{k), r{k) — r2(/c), where P2{k) = Pi{2i,i + I — k), 
q2{k) = {i + l - k){y - i - k), and r2{k) = {i + I + k){-2y -3i-l + 2j + k). So, here we 
have to find a polynomial S2{k) satisfying the recurrence 

Pi{2i,i + l-k) = {i + l-k){y-j-k)s2{k + l)-{i + l + k){-2y-3i-l + 2j + k)s2{k). (23) 
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Again, this is just routine for each specific instance of i and I, but we do not know an exphcit 
formula for S2{k) in general. Assuming that we have found a polynomial S2{k) satisfying 
(23), by (16) and (20) we have 



k=0 



r2(^ + / + l)g2(^ + / + l) , , , , r2(0)32(0) , 

t2{l + 1 + 1) T-T 12{0) 

P2{l + l + l) P2{0) 

{i + 



{i + l)\ 



i-2y -3i-l + 2j)n+i {y-j + l)n S2(0), (24) 



the last line being due to the fact that also t2{i + I + I) =0. 

In order to relate S2{k) to si{k), we make the following observation: We set S2{k) = 
S2{—k + 1), substitute this in the recurrence (23), then replace k by —k and change the sign 
on both sides of (23). Thus we obtain for S2{k) the recurrence 

-Pi{2i,i + l + k) = {i + l-k){-2y-3i-l + 2j-k)s2{k + l)-{i + l + k){y-j + k)s2{k). (25) 

This is almost the same recurrence as the recurrence (21) for si{k)\ It is only the term on 
the left hand side which is different! But, in fact, there is no difference: We claim that: 

Claim 1: We have P/(e, e + 2l- f) = -Pi{e, /). 

Claim 2: There exists a unique solution for the recurrence (21). 

Let us for the moment assume that these claims have been already established. Then, 
because of Claim 1, the recurrences (21) and (25) are indeed the same. Furthermore, thanks 
to Claim 2, there does exist a unique solution for the recurrence (21), and so also for (25). 
Hence, the solutions must be the same, i.e., si{k) = S2{k), which means si{k) = S2{1 — k). 
In particular, we have si(l) = ^2(0). A further fact, which follows immediately from Claim 1 
on replacing e by 2e and setting / = e + Z, is that P/(2e, e + I) = 0. Therefore, by setting 
A; = in (21), we obtain 

= (z + l){-2y -3i-l + 2j>i(l) -{t + l){y- j)si(O). 

From this equation, and the previous observation that si(l) = S2(0), we infer 

{-2y -3i-l + 2j) _ {-2y -31-1 + 2}) _ 

'■'"^ = = 

Substitution of this relation in (22) gives 

Z t / - fclr'"'' i + l-k) {-2y -3i-l + 2j-k + l)n+k {y + k - j + l)n-k 

k=0 ^' 

_ iz + -3i-l + 2j),+i {y-3 + 1)^^2(0). 



{i + l)\ 
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Comparison of this identity with (24) shows that indeed the sums on both sides of (14) are 
identical. This would prove (14). 

So it remains to settle Claims 1 and 2. 

We begin with Claim 1. By the definition (12) of Pi{e, /), we have 

21+1 

Pi{e, e + 2Z - /) = ^ (e)^ (-e -21 + f)2i+i-r 

r=0 
21+1 

= J2^r ie)r i-iy+' (e + r- f^l+l-r, 
r=0 

where the coefficients ar are given by (13). Next we use the Chu-Vandermonde summation 
(see e.g. [10, §5.1, (5.27)]) in the form 



XI \) (^)« iy)N-s = {x + y)N, 



s=0 

with N — 21 + 1 — r,x — e + r, and y — — /. Thus, 

21+1 2l+l-r 



PKe,e + 2/-/)= 5^a,(e),(-l)'^+i J] ( j (e + r), (-/)2^+i- 

21+1 m , , I 1 _ \ 

^-^ ^-^ \ m — r I 

^ — n 1 — n ^ / 



m=0 r=0 

Therefore, Claim 1 will follow immediately, if we are able to show that 



E(^^^^ M(-l)^a, = a^. (26) 
^-^ \ m — r I 



r=0 



This can be readily done by using generating functions. The definition (13) of the coefficients 
ar is equivalent to 

oo 

^arx"" = (a;^ +a; + l)'-^(2a;+ l)(a; + 2)(a;- 1). (27) 

r=0 

Let us denote the right hand side of this equation by A{x). Now we multiply both sides of 
(26) by a;"^, and we sum over all m = 0, 1, . . . We obtain 



m=Or=0 ^ ^ 



and after interchanging summations on the left hand side and summing the (now) inner sum 
over m by means of the binomial theorem. 
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It is trivial to verify this equation. Thus also the equivalent equation (26) must be true. 
Due to the preceding considerations, this completes the proof of Claim 1. 

Next we turn to Claim 2. We show that there is a unique polynomial si{k) of degree 21 
that satisfies the recurrence (21). (We leave it as an exercise that the "degree calculus" of 
the Gosper algorithm shows that if there is a solution to the recurrence (21) then it has to 
be a polynomial of degree at most 21.) So, let si{k) — X^^=o c(to)(/^ — ^ — Om- We substitute 
this into (21), then expand everything with respect to the basis {k — i — l)m, m = 0, 1, . . . 
(for the space of polynomials in k), and finally compare coefficients of {k — i — l)m on both 
sides of (21). This leads to the following system of equations for the coefficients c(m): 

a2i+i-m{2i)2i+i-m = {y + i - 1 - j + m - 1) c{m -l)-{2i + 2l-m){y + i + l-j-m) c{m), 

m = 0, 1,...,2Z+ 1, (28) 

where, by convention, we put c( — 1) = c(2/ + 1) = 0. For convenience, we set 

c(m) = {2i)2l-m (y + i-l- j)2l-m {y + i-l- j)m c{lT^)- 

By substituting this in (28), we obtain the simpler system of equations 

-— , ^ ■ , — = c(m-l)-c(m), m = 0,l,...,2/ + l. (29) 

[y + i-l- j)2i-m+i [y + ^ - I - J)m 

This is a system of 2Z + 2 equations for 2/ + 1 variables. (Recall the convention c(— 1) = 
c{2l + 1) = 0, which of course implies c(— 1) = c{2l + 1) = 0.) So, it is over deter mined. It 
is easy to see that this inhomogeneous system of linear equations has a (unique) solution if 
and only if the sum of the left hand sides of (29) over all m equals 0, i.e., if and only if 

21+1 

^ ^ (20) 

{y + i-l- j)2i-m+i {y + i-l- j)m 

This would follow immediately from the antisymmetry property = —a2i+i-rnj because 
then the m-th and (2/ + 1 — m)-th summand in the sum in (20) would cancel each other. 
Indeed, the substitution x 1/x in (27) yields am = —a2i+i-m- Therefore, the system of 
equations (28) has indeed a unique solution, which implies that there is a unique polynomial 
si{k) satisfying the recurrence (21), which is exactly the assertion of Claim 2. 

The proof that HiLil^ + 2y + 3i + l)n-i divides Dn{x, y) is now complete. 

Step 4- Dn{x,y) is a polynomial inx of maximal degree n{3n + l)/2, and the same is true 
for the maximal degree of Dn{x, y) in y. This is because each term in the defining expansion 
of the determinant Dn{x, y) has degree n{3n + l)/2 in x, and the same in y. Since the right 
hand side of (4), which by Steps 1-3 divides Dn{x,y) as a polynomial in x and y, also has 
degree n(3n + l)/2 in x, respectively y, Dn{x,y) and the right hand side of (4) differ only 
by a multiplicative constant. 

Step 5. The evaluation of the multiplicative constant. By the preceding steps we know 
that (5) holds. In particular, if we set y = 0, we have 

n 

det {{x + 2j-i+l)n+i{i-j + l)n-i) = C{n)T[{x + j + l)n+j. (31) 
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(The reader should be aware that the second term in the determinant Dn{x,y), as given 
by (4), vanishes for j/ = because of the presence of the factor {y — i — j + l)n+i-) The 
determinant on the left hand side of (31) is a lower triangular matrix, hence it equals the 
product of its diagonal entries, which is n^=i(^+J + l)n+j {n — jY-- Therefore C{n) is equal 

to n;=i(n-j)! = n;=i(j-i)!- 

This finishes the proof of (4) and thus of the Theorem. □ 

Acknowledgement. We are grateful to the referee for a simplification of our original proof of 
(26). 
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